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INTRODUCTION 

The control design problem for the class of future spacecraft referred to as 
large space structures (LSS) is by now well known [1-31. The issue is the re- 
duced order control of a very high-order, lightly damped system with uncertain 
system parameters, particularly in the high frequency modes. This paper pre- 
sents a design methodology which incorporates robustness considerations as 
part of the design process. Combining pertinent results from multivariable 
systems theory and optimal control and estimation, LOG eigenstructure assign- 
ment [4] and LQG frequency-shaping, [5-7] were used to improve singular value 
robustness measures in the presence of control and observation spillover. 

The design technique is summarized as follows. A low order LQG compensator is 
synthesized using the technique of recursive eigenstructure assignment to 
place closed-loop eigenvalues where desired. This design is evaluated for 
singular value performance margin and for singular value gain margin with 
respect to plant uncertainties (e.g., modeled dynamics). The compensator is 
then resynthesized using frequency-shaping concepts to improve the singular 
value robustness measures. The recursive eigenstructure assignment technique 
allows regulator close-loop eigenvalue placement at the desired locations for 
the plant and as required for frequency-shaping. Furthermore, the frequency- 
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shaped compensator eigenvalues can also be assigned, thus assuring LOG com- 
pensator stability, as well as estimator stability. 

This procedure using robust frequency-shaped compensation was applied to the 
design of the controller for a representative large space structure. Results 
are presented as singular value Bode i.!ots. Comparisons are made to a recent 
study9 utilizing the same large space structure model. 


LQG CONTROL DESIGN FOR LSS 

Control design plant modelling for LSS utilizes a high-order structural model, 
typically obtained by finite- element programs such as NASTRAN. The limita- 
tions of computer implementation require that the finite-element model be 
reduced to a design model. One approach is to truncate the high-order model 
into primary and residual modes, where the primary modes are to be used for 
control design. The modal truncation can be based on engineering judgement or 
on a selection criterion such as modal cost analysis [93. 

The system model has the form 

x p * Ap Xp + BpU 
• 

X R ” Ar xr + Bru (1) 

y - Cp Xp ♦ Cr xr 
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where x p are the primary modes and xr are the residual modes. An ob- 
server-based control design for v he primary modes then has the form 

X P » Ap x p + B p u + G(y - C p 5p) (2) 

A 

u = -KXp 

Using LQG design, the gains (K, G) are selected to minimize quadratic perform- 
ance indices. The terms Br u and Cr xr were identified by Balas [33 as control 
spillover and observation spillover respectively. These terms have the poten- 
tial fur interacting through the observer (2) to produce instability. 

LQG theory guarantees that the reduced-order closed loop system is stable with 

eigenvalues of (A p -B p K) and (A p -GC p ). However, no such guarantee holds for 
the compensator, 

u *= Hy (3) 

which has the eigenvalues of (A p -B p K-GC p ) . This fact can be fatal for LSS 
reduced-order control, unless measures are taken to ensure system robustness. 

ROBUSTNESS MEASURES FOR LSS 

For multivariable feedback systems the emerging singular vali *> robustness 
theory can be used to develop measures for stability and performance. Kosut, 
et al,8 applied this theory to the large space structure control design prob- 
lem, treating the residual dynamics as a perturbation. For a system with a 
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stable nominal feedback system (based on the reduced model) and stable pertur- 
bations (due to the residual dynamics), sufficient conditions for stability 
are obtained when the singular value stability measures exceed the maximum 
perturbation due to model uncertainty. Fig. 1 defines the terminology for a 
large space structure control system. For an additive perturbation, cig. 2a, 
the sufficient conditions for stability are 


SMi = ^ [I + H(jto) G c (ym)] > ^ [H(jo>) Gr (jto)] 
3M 2 = o [I + G c (jo>) H(ju)] > o [Gr (jto) H(jo,)] 


(4) 


where o (•) indicates the maximum singular value and o_ (•) indicates the mini- 
mum singular value. (Singular values of the complex matrix A are the positive 
square roots of the eigenvalues of A*A, where (•)* indie n os conjugate 
tran-spose.) If G c (s) is minimum phase and invertible, a multiplicative per- 
turbation can be formed, Fig. 2b, and the sufficient conditions for stability 
are then 


SM^ = _o [I + (HG C )~1] > o [Gq-^Gr] 


(5) 


S.1 ? = £ [I + (G C H) -1 ] > o [GrGc' 1 ] 
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where the ju arguments have been supressed. Good performance within the oper- 
ating frequency region (i.e., the "control bandwidth") is provided when the 
performance measure 

PM « o [I ♦ G c h] (5) 

is large. The stability measures (i») are generalizations of Nyquist polar 
plot analysis; the measures (5) are generalizations of Nyquist inverse polar 
plot analysis. The need for large performance measure (6) is a generalization 
of the desirability of large loop gains. 

ROBUST COMPENSATION DESIGN 

The stability and performance measures presented above require stac, lity of 
the nominal feedback system. In a previous work £*1], the authors presented a 
recursive design procedure which assigns the closed-loop eigenstructure in 
linear quadratic regulators. At each stage, the required solution for the 
steady state hiccati matrix which shifts a pole or po A e pair to specified 
values is obtained. For pole pair placement, a free parameter in the solution 
permits selection of closed-loop eigenvectors. This design procedure is sum- 
marized in Appendix 1. 

Using duality, the procedure also applies to estimator design. By extension, 
the procedure can be used to design stable compensators by considering the 
closed-loop regulator dynamics matrix (A-BK) as the open-loop system and pick- 
ing the estimator gain to place the compensator eigenstructure of (A-BK-GC). 
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Compensator robustness can be enhanced through the use of frequency-shaped 
control and estimation [5,6]. In frequency-shaped estimation, a fre- 
quency-domain perf ormance index is considered, 


1 00 

J = E — / [w'Q(ju)w + v' R(jw)v] du> (7) 

2n -* 

where w is the disturbance and v is the sensor noise. Sensor r.oise fre- 
quency-shaping is re'll ized by treating v as an auto correlates noise source of 
the form 


v(j w) 


= R 



u) V(ju) 


(3) 


where v'(ju) is a white noise process. In the approach used here, Q(jui) is 
determined by pole placement, equivalent to injecting fictitious process 
noisj. Rl/2(j u) ) must be proper (not strictly proper) to maintain sensor noise 
weighting over the entire spectrum. Then define a pseudo-measurement 

z' - R-1/2 (j u ) z - R-1/2(J U ) Cx(Jid) ♦ v' (ju) (9) 

R" 1// 2 (ju) can be realized in state space as 

J 

Xy ■= Ay Xy + ByCX 

y * Cy Xy + DyCX (10) 

2’ * CyXy + DyCx + V* 
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This dynamic model is appended to the system dynamics to form the frequency- 
shaped estimator, 

x * Ax + G(z’ - C V X V - D v Cx) + B c u 

( 11 ) 

4 

Xy * AyXy + ByCX + Gy (Z* “ CyXy — DyCx) 


where z' is obtained from (10). The gains G and G v t jn be kicked to plav'. the 
eigenvalues of (11) at those of the frequency-shaping filter (10) and the oth- 
ers as required for performance. A dual result can be used to develop fre- 
quency-shaped gains for the regulator. 

Because frequency-shaping ada... states to the compensator, an efficient choice 
of the loops to be shaped Is desirable. Kim [7] has developed a procedure for 
loop selection based on the singular vectors of the return ratio matrices G C H 
or HG C . Me conjectured that an input vector y in the direction of qj , the 
singular vector corresponding to "oik) will get the largest amplification by A. 
Similarly, a vector in the direction of q n , the singular vector corresponding 
to _o(A) will get smallest amplification. Therefore, if the component of y In 
the direction which is closest to qi , is reduced by a filter before it enters 
A, ~(A) is effectively reduced. a( A) increased by increasing the compos 
of y closet tc q n before it enters A. It > -n be shown that frequency-shaping 
introduces transmission zeros into the compensator transfer function. 
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DESIGN METHODOLOGY 


The discussion which has been presented above suggests the following design 
methodology : 

. Compensator design for performance of the reduced order system. 

2. Eva* .ation of the stability margins (4,5) against the perturbation due to 
the residual dynamics. 

3. Selection of frequency-shaping filters to enhance stability robustness. 

4. Synthesis c: . equency-shaped compensator to incorporate performance and 

stability margins. 

The recursive eiger.structure design algorithm can bo used for the designs. 


EXAMPLE. 


Tne design methodology was applied to a control design for the ACOSS-1 model, 
albj u'.ed in the comparison study [9]. The model is illustrated and the 
state-space data are listed in Appendix 2. As in the comparison study *'he 
first eight structural modes were retained. A regulator was designed with 
ciosed-loop poles at 20$ damping; a compensator was designed with poles at 
ritical damping. Fig. 3 illustrates stability measure (5) for the loop 
broKen at the ouunut. Performance s adequate at low frequencies but stabil- 
ity robustne:- is inadequate above 1 Hz. 
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To improve stability robustness, frequency-shapea estimation was incorporated 
in all three output loops using second-order low-pass filters. Fig. 4 illus- 
trates the recovery of stability robustness while still retaining good low 
frequency performance, Fig. 5. 


DISCUSSION 


In the comparative study by Kosut, et al [8], both LOG modal control and a 
frequency-shaped control were investigated (along with others). LOG control 
was found to have poor performance as well as poor stability robustness. Fre- 
quency-shaped control was found to have adequate stability robustness, but 
poor low frequency performance. 

The methodology presented here addresses both of these issues. Performance is 
achieved by pole placement design of the compensator, achieving good loop 
gains at low frequency. Stability robustness is achieved by adding fre- 
quency-shaping without sacrificing low frequency performance, since the gain 
of the frequency-shaping filters is one at low frequencies. 


CONCLUSIONS 


A design meth dology for control systems for large space structures has been 
proposed which incorporates both performance and stability robustness concerns 
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as an integral part of the design process. Performance was achieved by plac- 
ing the poles of the compensator. Stability robustness was achieved by fre- 
quency-shaping the compensator to satisfy a frequency domain stability 
robustness test. 

An example was pr isented which applied the methodolgy to a system with the 
loop broken at the output. A full design study would also require examination 
of the system with the loop broken at the input, using regulator fre- 
quency-shaping to enhance robustness. 
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APPENDIX 1 


Recursive Eigenstructure Design 

The steady-state optimal control law for the linear, time-invariant, 
controllable system: 

x = Ax + Bu (A.1) 

which minimizes the quadratic performance index, 

J = 1/2 lxT q x + U T Ru] dt (A. 2) 

is linear state feedback 

u * Kx = - R'l bT sx (A. 3) 

where S is the solution of the steady-state Riccati equation, 

-SA - AT S + SBR-1 BTS - Q - 0 (A. U) 

In this appendix we summarize an interactive design technique which solves 
(A. 4) to provide specified eigenvalues of the closed-loop system dynamics 
matrix A+BK and which also permits some freedom in selecting closed-loop 
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eigenvectors. The method is reported elsewhere [4] in detail. It extends the 
procedure of Solheim [10] in which, for fixed R, the elements of Q providing 
the required pole placement are calculated directly. 

The design technique is recursive; at each stage, the system dynamics matrix A 
in ( A. 1 ) incorporates previous state feedback. We then implement the follow- 
ing eigenstructure calculation: 

X"1[A - HS]X = A (A. 5) 

where A = T -1 AT is block diagonal, T is the real eigenvector matrix of A, and 
H = T~1 br _1 B t T~ t is symmetric and positive semi-definite. A is identical to A 
except for a block of shifted poles. X is the transformation from open-loop 
eigenvectors to closed-loop eigenvectors; it is defined as the "stage" 
eigenvector matrix. ? is the Riccati matrix in the open-loop diagonalized 
coordinate system; S is chosen to shift a single pole or a pair of poles. The 
corresponding gain matrices, K, de ermined for each stage are subsequently 
added to obtain a final gain which achieves the same closed-loop pole loca- 
tions. 

To provide the required pole shift, the only non-zero elements of S correspond 
to the entries of A which are to be shifted. With this ohoice of S, the char- 
acteristic equation factors into the product of terms for the unshifted poles 
and a term for the desired shifted poles. Thu~, 

f si - a! - D(s) n (s - Ai) (A. 6) 

iel 
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where I is the index set for the unshifted poles, and D(s) contains explicit 
elements of S, H, and the block of A which is to be shifted. Matching the 
coefficients of powers of s in D(s) to the equivalent terms in the closed-loop 
characteristic equation provides a set of equations in the required elements 
of S. For the single pole shift A jj - * to ic, the only non-zero element of S 
satisfies 


For double pole placement it can be shown that the three required elements of 
lie on the intersection of two quadric surface in a mathematical space hav- 
ing the three S elements as coordinates. (It can also shown that a direct 
solution for Q has a similar geometric interpretation.) If the corresponding 
submatrix of H is positive definite, the surfaces are a plane and a hyper- 
boloid of one or two sheets; the inters:;' -Ion, if it exists, is always an el- 
lipse. If the relevant submatrix of L is singular, the surfaces are planes, 
and the intersection is a line. The different points comprising the solution 
all provide the desired eigenvalue placement, but with different eigenvectors. 

In ref. ^ a solution for S is presented which takes advantage of the quadric 
surface geometry to define a free parameter that allots design freedom in the 
choice of closed-loop eigenvectors. The solution for the stage eigenvector X 
partitions into two sets of equations. The first is a homogeneous Lyapunov 
equation for the submatrix corresponding to the shifted pole block in A. For 
> role pair shift, the submatrix is 2x2. Hence, depending upon the nature of 
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the closed-loop poles (real or complex), one or two elements of the submatrix 
may be chosen arbitrarily; the remaining elements then depend :-n the choice of 
elements of S. The other equation is a non-homogeneous Lyapunov equation in 
the remaining elements of the columns of X containing the 2x2 submatrix; its 
solution depends upon the 2x2 submatrix, the elements of S and A, and certain 
elements of H. 

The closed-loop system eigenvector matrix is then Tql - TX. The solution of X 
depends upon S, which varies with the choice of the free parameter. There- 
fore, by recursively shifting pole pairs, design freedom exists to select 
closed-loop eigenvectors while providing required pole placements. 

The procedure outlined above lends itself to a recursive procedure for practi- 
cal multivariable regulator design. The steps in the procedure are as fol- 
lows: 

1. System (A.1) is placed in modal form. 

2. The designer selects the control weighting matrix R, then H is calculated. 

3. The designer selects a real pole or pair of poles to be shifted and their 
desired location; • a pair, he also selects the free parameter which 
determines the closed-loop eigenvectors. 

4. The stage gain is calculated and the closed-loop system is placed in modal 
form. 
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5. Steps 3 and 4 are repeated for other poles until the designer is satis- 
fied. 

6. The total system gain is obtained by adding the stage gains. 

Clearly by duality, the same process can be applied to estimator design, per 
mitting the development of multivariable compensators. 


78 



APPENDIX 2. 


The ACOSS-1 flexible spacecraft model was developed by the Charles Stark 
Draper Laboratory . 10 it is representative of many radar and optical control 
problems, but is small enough to be tenable for research studies. The struc- 
ture is a tetrahedral truss supported by three right-angle bipeds. The truss 
members are flexible in the axial direction only. The model has 12 modes; for 
control design, only eight are assumed to be -known. 
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Residual Dynamics 
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